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ABSTRACT 

M 87 is the only known non blazar radio galaxy to emit very high energy (VHE) gamma-rays. During a 
monitoring program of M 87, a rapid flare in VHE gamma-rays was detected by the MAGIC telescope in early 
2008. The flux was found to be variable above 350 GeV on a timescale as short as 1 day at a significance level 
of 5.6a. The highest measured flux reached 15% of the Crab Nebula flux. We observed several substantial 
changes of the flux level during the 13 day observing period. The flux at lower energies (150 - 350 GeV), 
instead, is compatible with being constant. The energy spectrum can be described by a power law with a 
photon index of 2.30 ±0.1 l sta t±0.20 S y St . The observed day-scale flux variability at VHE prefers the M 87 core 
as source of the emission and implies that either the emission region is very compact (just a few Schwarzschild 
radii) or the Doppler factor of the emitting blob is rather large in the case of a non expanding emission region. 

Subject headings: gamma rays: observations — galaxies: active — galaxies: individual (M 87, NGC 4486) 
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1. INTRODUCTION 

M 87 is a giant elliptical radio galax y (RG) of Fanaroff- 
Riley-I-type (FR I, iFanaroff & Rilevll 19741) in the Virgo Clus- 
ter at a distance of 16Mpc dMacri et al.lll999l) . It is powered 
by a supermassiye blac k hole (BH) of (3.2 ± 0.9) x 10 9 M Q 
(iMacchetto et al.lll997l) . The M 87 jet was the first-ever ob- 
served (ICurtislll9 181) 7 and due to the proximity of M 87, its 
morphological substructures can be resolved and a unique 
view of its innermost regions is possib le. The jet, originat- 
ing from the RG core, extends to 20" dMarshall et al.1 120021: 
equivalent to a 2 kpc projected linear distance). Several com- 
pact regions ("knots") along its axis are resolved in the radio, 
optical, and X-ray regimes. These knots have similar mor- 
phologies in all wave bands, although the X-ray knots appear 
to be t ens of pc closer to the core than the optical and radio 
knots (Wil son & Yandl2002l) . The variable brightness of the 
knots may be due to several shock fronts in the jet, being re- 
sponsible for particle acceleration and non-thermal emission. 
Superlu minal motion of the knots has been observed in the 
optical (lBirettal["l996l) and radio dForman et al.l l2007h wave 
bands, constraining the viewing angle of the jet to < 43° ±4°. 

The innermost resolved bright knot HST-1 is located at 
0.82" (64 pc) from the core and is the most prominent feature 
of the jet. HST-1 has shown many flares exceeding the lu- 
minosity of the M 87 core emission. Its X-ray brightness has 
increased by a factor > 50 from 2000 to 2005 (lHarrisll2006l) . 
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A correlation between radio, optical, and X-ray luminosity 
points to a common origin of the emission. The measured 
superluminal motion in HST-1 is higher than in other knots, 
suggesting a viewing angle of < 19° for this part of the jet. 
The core itself is variable, too, and also shows a correlation 
betwee n the emission levels from radio frequencies through 
X-rays dPerlman et al.ll2003l) . 

M 87 was not detected by EGRET. The first hint of very 
high ener gy (VHE; > 25 GeV) 7-ray emission was re- 
ported by Aharonia net al.l (120031). and l ater c onfirmed by 
lAharonian et all d2006l) and lAcciari et all (120081) . The emis- 
sion is vari able on a timescale of y ears. The reported « 2-day 
variability ( Aharoni an et al.ll2006l) narrows down the size of 
the emission region to be on the order of the light-crossing 
time of the central BH. With its expected low accretion rate, 
the M 87 core radiation is not strong enough to attenuate sig- 
nificantly Te V 7-rays even at 5 Schwarzs child radii (Rs) away 
from the BH dNeronov & Aharonianl2007h . All this implies a 
production region in the immediate vicinity of the M 87 core. 
Durin g later observations , no significant flux variation was 
found ( lAcciari et al.ll2008l) . An X-ray-VHE 7-ray correlation 
is expected in most emission models, but was not un ambigu- 
ously found so far. Whereas Ahar onian et al.l (12006b claim a 
hint of a correlation betwee n the soft (0.3-10 ke V) X-rays at 
HST-1 and the VHE -7-ravs. lAcciari et al l d2008l) find a year- 
by-year correlation between the 2-10 keV X-ray flux of the 
M 87 core and the VHE 7-ray emission instead, but do not ob- 
serve a correlation between the two energy bands on shorter 
timescales. 

The radio to X-ray emission of the jet is due to non-thermal 
synchrotron radiation of relativistic electrons in the jet. The 
observed knots and flares in M87 point to a complicated 
morphology with several shock fronts producing these elec- 
trons. While all 23 currently known extragalactic VHE 7- 
ray emitters 30 are blazars, M87 i s assumed to be a blazar 
not aligned to our line of sight dTsvetanov et al.l fl998l) . If 
the observed VHE emission from M 87 is associated with the 
innermost part of its jet, then blazar emission models may 
hold. In blazars without prominent disk or broad-line fea- 
tures, the VHE emission is explained by inverse Compton 
processes involving the synchrotron photons and their par- 
ent electron population (synchrotron self-Compton models, 
e.g. lMaraschi et al1 fl992). Alternatively, in hadronic models, 
interact ions of a highly relativistic jet outflow with ambient 
matter dDar & LaoiHl997t iBeall & Bednareklll999l) . proton- 
induced cascades (Mannheim 119931) , or synchrotron proton 
radiation (Mucke & Protheroel kOQll: lAharonianl l2000h may 
produce VHE 7-rays. In such a scenario, M 87 might also 
accou nt for parts of the o bserved ultra-high-energy cosmic 
rays dProtheroe et alJ 120031) . It should be noted, that for 
M87 the location of the VHE emission is still uncertain. 
Specific emission models for high-en e rgy processes close to 
the core (iGeorganopoulos et al J 120051: iGhisellini et al.l 120051: 
iLenain et al.l 120081: iTavecchio & Ghisellinil 120081). in th e 
large-scale jet dStawarz et alJ 12003b iHonda & Hondall2007l) 
or in the vicinity of a B H dNeronov & Aharo nian 2007; 
iRieger & Aharo nian 2008) have been developed. 

The MAGIC collaboration performed monitoring observa- 
tions of M87 starting from 2008 January, sharing the task 
with the VERITAS and H.E.S.S. experiments. Here we re- 
port on MAGIC results from a subset of these data, revealing 

30 See, e.g., http://www.mppmu.mpg.de/~rwagner/sources/ for an up-to- 
date list 



a variability timescale of M 87 of 1 day. The energy spectrum 
and light curves are discussed. 

2. OBSERVATIONS AND DATA ANALYSIS 

The MAGIC telescope is located on the Canary Island of La 
Palma (2200 m above sea level, 28°45'N, 17°54'W). MAGIC 
is a stand-alone imaging air Cerenkov telescope (IACT) with 
a 17 m-diameter tessellated reflector dish. MAGIC has a low 
energy threshold of 50 - 60 GeV (trigger threshold at small 
zenith angles). The accessible energy range extends up to tens 
of TeV with a typical energy resolut ion of 20%-30%, de pend- 
ing on the zenith angle and energy d Albert et al.ll2008a) . 

The data set comprises observations from 2008 January 30 
to 200 8 February 11. Th ese were performed in the wobble 
mode (Da um et al.l fl997) for 26.7 hr. The zenith angle of 
the observations ranges from 16° to 35°. After removing 
runs with unusually low trigger rates, mostly caused by bad 
weather conditions, the effective observing time amounts to 
22.8 hr. 

The data were a nalyzed using the M AGIC standard calibra- 
tion and analysis dAlbert e t al. 2008 j). The analysis is based 
on image parameters jlHillasI IT"985b iTescaro et alJl2007l) and 
the random forest (RF, Alber t et al . 2008b) method, which are 
used to define the so-called hadronness of each event. The 
cut in hadronness for 7/hadron separation was optimized on 
a contemporaneous data set of the Crab Nebula. After this 
cut the distribution of the angle ALPHA, which is the angle 
between the main image axis and the line between center of 
gravity of the image and the source position in the camera, is 
used to determine the signal in the ON-source region. Three 
background (OFF) sky regions are chosen symmetrically to 
the ON-source region with respect to the camera center. The 
final cut \ALPHA\ < 5° (Fig. [T]) was also optimized on the 
Crab Nebula data to determine the number of excess events 
and the significance of the signal. 

The energies of the 7-ray candidates were also estimated 
using the RF method. To derive a differential energy spec- 
trum, we applied looser cuts than those in Fig. Q] to retain 
a higher number of 7-ray candidates and to lower the effec- 
tive analysis threshold down to 150 GeV. Looser cuts also re- 
duce systematic uncertainties between data and Monte-Carlo 
events, which is important for the estimation of the effective 
collection areas. The derived spectrum was unfolded to cor- 
rect fo r the effects of the limited energy resolution of the de- 
tector (Albert et al. 2007). Finally, the spectrum and the light 
curves were corrected for trigger inefficiencies due to higher 
discriminator thr esholds during par tial moon light and twi- 
light conditions dAlbert et alj|2008cl) . These corrections are 
on the order of 0%-20%. The data were also analyzed with 
an independent analysis yielding, within statistical errors, the 
same results. 

3. RESULTS 

The \ALPHA \ distribution after so-called tight cuts is shown 
in Fig. [T] The applied cuts are SIZE > 450photoelectrons and 
hadronness h < 0.02. After the final \ALPHA\ cut (resulting 
in an overall cut efficiency of 37%), the total signal of 241 
events over 349 normalized background events corre sponds to 
an exc ess with a significance of 9.9ct along eq. 17 in lLi & Mai 
d!983l) . The highest flux was observed on 2008 February 1 at 
a significance of 8. Oct. 

In searching for time variability, the data set was divided 
into nine subsets, one per observing night. In Fig. [2] we 
show both the light curve above the analysis threshold (150- 
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M87, January-February 2008 
Obs.Time: 22.8h 
Significance: S = 9.9 a 
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FIG. 1. — \ALPHA\ distribution for the overall data sample. The back- 
ground (gray histogram) is estimated using three OFF regions arranged sym- 
metrically to the ON-source region with respect to the camera center. A 7-ray 
excess with a significance of 9.9 standard deviations is obtained. 
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FIG. 3. — The differential energy spectrum of M 87 for the total data sample. 
The horizontal error bars represent width of the energy bins. The best-fit 
function, assuming a power la w, is given by the solid curve. The Crab Nebula 
spectrum ( Albert et al. 2008a) is given by the dashed curve for reference. 
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FIG. 2. — The night-by-night light curve for M87 as measured from 2008 
January 30 (MJD 54495) to 2008 February 1 1 (MJD 54507). The upper 
panel shows the flux in the energy bin 150 — 350 GeV, being consistent with 
a constant emission. The lower panel shows the integral flux above 350 GeV; 
flux variations are apparent on variability timescales down to 1 day. The inlay 
of the lower panel shows the light curve above 350 GeV in a 40 min time 
binning for the night with the highest flux (2008 February 1). The vertical 
arrows represent flux upper limits (95% c.l.) for the nights with negative 
excesses. 



350 GeV) and in the energy range at which MAGIC has the 
highest sensitivity for a variability search (>350GeV). The 
low-energy range shows no significant variability with a xt 
of 12.6/8 (probability of P = 0.13) for a constant fit. In- 
stead, in the energy range above 350 GeV clear variability is 
found. A fit by a constant has a xl of 47.8/8 corresponding to 
P=l.l ■ 10~ 7 . The correlation coefficient between the two en- 
ergy bins is r = — 0.25^33 (1-cr errors), suggesting that there is 
no significant correlation, but we note rather large error bars in 
the low energy flux bin. We also investigated a night-by-night 
variability. There are five pairs of observations on consecutive 
nights in the total data set. We calculated individual probabil- 
ities St for these pairs to have the same flux level and the cor- 
responding signi ficances. We then co mputed a combined sig- 
nificance S com b (Bit vukov et alj|2006» : S CO mb = (X^/) /s/n, 



with n = 5. We interpret the resulting S com b = 5.6 a as a proof 
that the flux varies on timescales of 1 day or below. Note that 
the 1 day variability is claimed from this combined analysis 
rather than from the 2008 February 1 flare alone. We find our 
statistics not sufficient enough to determine the flare shape. 
Given the number of the observed changes in the flux level, 
the data belong to a complex of two, if not three, sub-flares. 

We also looked for shorter time variability, but in none of 
the observation nights there is a significant flux variation in 
the two energy bands. A typical example in a 40 minute bin- 
ning is shown in the inset in Fig.|2]for 2008 February 1. 

The averaged differential energy spectrum of M 87 (Fig. [3) 
extends from ^100 GeV to ~ 10 TeV and can be well ap- 
proximated by a power law: 
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The errors are statistical only. We estimate an 11% sys- 
tematic uncerta inty in the normal ization and 0.20 for the 
spectral index (Alb ert et alj l2008ab . The measured values 
are in good agr eement with the H.E .S.S. (spectral index 
T = -2.2 ± 0. 15. lAharonian et al.l l2006l) and VERITAS (T = 
-2.31 ±0.17. lAcciari et al.E 008) results. The observed spec- 
trum is not significantly affected by the evolving extragalac- 
tic background light (EBL) due to the proximity of M 87 
dNeronov & Aha ronian 2007). To investigate a possible hard- 
ening of the spectrum with increasing absolute flux level, we 
divided the data sample into high and low state subsamples. 
The high sample comprises the two nights with the highest 
flux above 350 GeV (February 1 and 8), while the low state 
comprises the nights of lower-flux data (January 30, February 
2, 4, and 1 1). Both the high and low state spectra (Fig. |4j can 
be well described by a power law: 
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(4.81 ±0.82), r hi s h = (-2.21 ±0.18) and /^ ow 



with /„ 

(2.58 ±0.92), T low = (-2.60 ±0.30) for the high and low 
states, respectively. There is a marginal hardening of the spec- 
tral index with the higher flux on the level of 1 - 2 standard de- 
viations, depending on the way the significance is calculated. 
This hardening is not significant, which can be a consequence 
of the fact that the two flux levels (states) differ by less than a 
factor of 2. 
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FIG. 4. — Differential energy spectra of M87 divided into high (filled 
circles) and low (open circles) states. See text for the details. The best-fit 
functions, assuming power laws, are given by the black solid and red dashed- 
dotted curves, respectively. 



4. DISCUSSION 

M87 is the only non blazar radio galaxy known to emit 
VHE 7-rays and one of the best-studied extragalactic black- 
hole systems. To enable long-term studies and assess the 
variability timescales of M87, the H.E.S.S., VERITAS, and 
MAGIC collaborations established a regular, shared monitor- 
ing of M 87 and agreed on mutual alerts in case of a significant 
detection. During the MAGIC observations, a strong signal of 
8(7 significance was found on 2008 February 1, triggering the 
other IACTs as well as Swift observations. For the first time, 
we assessed the energy spectrum below 250 GeV, where our 
observations can be well described by a power law that shows 
no hint of any flattening. 

Our analysis revealed a variable (significance: 5.6 a) night- 
to-night 7-ray flux above 350 GeV, while no variability 
was found in the 150-350GeV range. We confirm the 
E > 730 GeV short-ti me variability of M87 reported by 
lAharonian et al.l d2006l) . The observed variability timescale is 
on the order of or even below 1 day, restricting the emission 
region to a size of R < AtcS = 2.6 x 10 I5 cm = 2.6SR S . The 
Doppler factor S is only relevant for an emission region not 
expanding while traversing the jet. In case of an expanding- 
jet hypothesis, the initial radius of the expanding shell is given 
by R* < cAt. The emission can occur very close to the BH, 



provided that the ambient photon density is low enough as to 
allow the propagation of VHE 7 rays. Otherwise the emission 
region must be located farther away from the BH. In the latter 
case, the variability constraints can be met only if the emitting 
plasma does not substantially expand while traversing the jet, 
or if it moves with 5 > 10. 

There exists no lower limit on the size of HST-1, and thus 
the flux variability cannot dismiss HST-1 as possible origin 
of the TeV flux. During the MAGIC observations, however, 
HST-1 was at a historically low X-ray flux level, whereas at 
the same time the luminosity of the M 87 core was at a histor- 
ical maximum (D. Harris 2008, private communication). This 
strongly supports the core of M 87 as the VHE 7-ray emission 
region. 

Our data alone cannot put strong constraints on VHE 7-ray 
emission models. The relatively hard VHE spectrum found 
for M87 (r ps -2.3) is not unique among the extragalactic 
VHE 7-ray sources if one considers intrinsic spectra, i.e. EBL 
corrected. Also, we did not measure a high-energy spectral 
cut-off. The found marginal spectral hardening may be inter- 
preted as a similarity to other blazars detected at VHE, where 
such hardening has often been observed. 

Our results show that a dense TeV monitoring, as exercised 
by ground-based IACTs, has revealed highly interesting rapid 
flares in M 87. This fastest variability observed so far in TeV 
7-rays from M 87 restricts the size of the 7-emission region to 
the order of R$ of the central BH of M 87 and suggests the core 
of M87 rather than HST-1 as the origin of the TeV 7-rays. 
Results from the entire monitoring campaign, comprising data 
from other IACTs, will appear in a separate paper. 
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supported by ETH Research Grant TH 34/043, by the Pol- 
ish MNiSzW Grant N N203 390834, and by the YIP of the 
Helmholtz Gemeinschaft. We thank Dan Harris for providing 
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